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Abstract. The so-called solar cycle is generally characterized by the quasi-periodic oscillatory 
evolution of the photospheric spots number. This quasi-periodic pattern has always been an 
intriguing question. Several physical models were proposed to explain this evolution and many 
mathematical data analysis were employed to determine the principal frequencies noticeable in 
the measured data. Both approaches try to predict the future evolution of the solar activity and 
to understand the physical phenomena producing these cycles. Here we present the analysis of 
the sunspots number evolution using the time-delay approach. Our results show than the solar 
cycle can also be characterized by this behavior implying the influence of the past evolution over 
the present one, suggesting an histeresis mechanism, linked probably with magnetic activity. 
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1. Introduction 

Solar activity can be seen through the evolution of sunspots number in quasi-periodic 
oscillatory series with periods going from 8 to 15 years and with a mean period of 11 years. 
Due to the change of magnetic field polarity in solar hemispheres alternatively each cycle, 
the period is rather 22 years. The quasi-periodic evolution of many activity phenomena 
is still a unsolved key problem in solar physics (along with, e.g., heating of the solar 
chromosphere and corona, and solar flares). An important issue of this understanding is 
due to the influence of solar activity over the terrestrial climate (Archibald 2006). 

Many physical models were proposed to understand the basic mechanisms (e.g., Benev- 
olenskaya 1998) producing solar/stellar activity and its quasi-cyclic evolution. From 
sunspots time series, mathematical approaches have highlighted several other hidden 
periods other that the well-apparent 11-years period. These works contributed to pro- 
vide some observed parameters to constraint theoretical models, and to predict the future 
evolution of the sunspots number (e.g., Clilverd et al. 2003; Sello 2003). 

In this work we follow another mathematical approach, investigating the influence of 
the past evolution of the sunspots number over the present one. A natural and simple 
way to carry out this analysis is to assume that the susnpots number have a temporal 
delay behaviour. Thus, the relationship between present and past values becomes non- 
linear and could prove that the past amplitude influence the present one. As for the 
abovementioned methods, this one also allows some predictions about future solar cycles 
by using intermediate parameters to characterize the general past evolution. 

2. The solar cycle as a temporal delay phenomena 

We consider a temporal delay behavior (Eq. 12. ip to rely the present evolution of a 
phenomena with some particular events in its past evolution. The variation of N in time 
is not only related with its current value N{t) but also with its past values N{t — T), 
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Figure 1. Left panel: Correlation 
factor obtained by chi-square mini- 
mization versus the temporal delay. 
Right panel: Sunspots number ratio 
versus sunspots number with a 8-years 
delay. The blue line shows a linear fit. 




where T is our temporal delay. The simplest way is to assume only one past temporal 
influence and some proportionality between the variables (e.g., Murray 1993), as follows: 



dN{t) 



a N{t){l-b N{t~T)} 



AN,, 



(2.1) 



dt ' ' ' V /J 

Because we use yearly values of the relative sunspots number, given by SDIC (http:// 
www.sidc.be/sunspot-data/). At is equal to 1. The present and past values of N are 
already known and we assumed some values for the parameter T {T — 0,1,2, . . .; Fig. [H 
left panel). So we can determine constants a and b fEa. l2.2p from the correlation between 
the sunspots number ratio and the past sunspots number (Fig. [T] right panel). 
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3. Results and conclusions 

The T values which minimize the correlation between the sunspots number ratio and 
the past spots number, are 7 or 8 years (Fig. [T|). To test the accuracy of our method, 
we applied it on the solar cycle 23 (Table [T] and Fig. [21 left panel). Our predictions are 
consistent with the observations, except for the unusually long period of low activity at 
the end of this cycle. We also characterized the next solar cycle 24 (Table [2] and Fig. [2]) 
using our modelling of the solar cycle 23 as well as the spots number observed until 
2006.5 and 2007.5. The next maximum sunspots number would take place between 2011 
and 2012 and should be close to 60 (Fig. [21 right panel) as for the solar cycle 14. Using 
the solar minimum occurred in 2008.5 {Rmin = 2.9 from SDIC), the predicted maximum 
is consistent with the value from the linear least-square fit on the both axis (red line on 
Fig. [3]) and is a few smaller than the range given by Brajsa et al. (2009). Moreover the 
next solar minimum should not occur before 2019 or even 2020. 



Table 1. Observed and predicted values for the cycle 23. Results showed in Fig. [2] (left panel). 



Observed parameters 


Epoch of solar 


Maximum sunspot 


Epoch of the end 


of solar minimum 


maximum 


number 


of cycle 


Epoch Number 


Observed Predicted 


Observed Predicted 


Observed Predicted 


1996.4 8.0 


2000.3 2001.5 


120.8 106.4 


2009.? 2006.5 



These preliminary results are encouraging because we find a similar delay as that 
observed between the geomagnetic activity and solar cycle peaks (Hathaway & Wilson 
2006) linked probably with magnetic activity by some kind of histeresis mechanism. 
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Figure 2. The observed {plus symbols and dotted line) and selected {asterisks and blue line) yearly values 
of the relative sunspots number. Our predictions using the temporal delay method are also marked {diamonds 
and red line). Left panel: Prediction for the solar cycle 23 and 24. Middle and Right panels: Prediction of the 
solar cycle 24 using the sunspots number observed until 2006.5 and 2007.5 as input. 



Table 2. Predicted values for the solar cycle 24 using various data set in input. 



Panel of 


Input parameters of 


Predicted parameters 


Predicted parameters of 


Fig.E] 


solar 


minimum 


of solar 


maximum 


next solai 


minimum 




Epoch 


Number 


Epoch 


Number 


Epoch 


Number 


Loft 


2006.5 


19.6 


2011.5 


84.1 


2016.5 


23.0 


Middle 


2006.5 


15.2 


2012.5 


137.3 


2018.5 


28.8 


Right 


2007.7 


7.5 


2012.5 


57.1 


2020.5 


4.2 



Figure 3. Maximum (Rmax) versus minimum (Rmin) sunspots num- 
ber of a given solar cycle (dots). The cycle number is also marked. Blue, 
green and red diamonds show the locus of our predicted maxima using 
our modelling of the solar cycle 23 as well as the spots number observed 
until 2006.5 and 2007.5, respectively (Tablc[2}. The predicted range of 
Brajsa et al. (2009) is plotted as well [vertical line). The blue and red 
lines show linear fits using data of both axis and the linear least-square 
fit like one given by Brajsa et al. (2009, Fig. 6), respectively. 




At present, we can successfully reproduce some previous solar cycles (e.g., epoch and 
solar maximum). We plan to improve our predictive method including the influence of 
delays for which there is a good correlation. This more detailed description could allow 
to solve the overestimation of the next solar minimum and to reproduce the asymmetries 
observed during the previous cycles. Our final aim is to obtain a reliable prediction of the 
whole solar activity by identifying the fundamental lower and high activity precursors of 
its present and future evolution. All these improvements are still needed to be able to 
predict future low solar activity periods (e.g., Maunder/Dalton Minimum). 
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